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A Universe filled with a homogeneous scalar field exhibits ^Cosmological hystere- 



o 

^S| ' sis\ Cosmological hysteresis is caused by the asymmetry in the equation of state 

^ ' during expansion and contraction. This asymmetry results in the formation of a hys- 

^f) , teresis loop: <f pdV, whose value can be non- vanishing during each oscillatory cycle. 

For flat potentials, a negative value of f pdV leads to the increase in amplitude of 
consecutive cycles and to a universe with older and larger successive cycles. Such 



i_j. a universe appears to possess an arrow of time even though entropy production is 

cn ! absent and all of the equations respect time-reversal symmetry ! Cosmological hys- 
> 

•^ ' teresis appears to be widespread and exists for a large class of scalar field potentials 



and mechanisms for making the universe bounce. For steep potentials, the value of 



cn . S pdV can be positive as well as negative. The expansion factor in this case displays 

CN I quasi-periodic behaviour in which successive cycles can be both larger as well as 

smaller than previous ones. This quasi-regular pattern resembles the phenomenon 
/\ . of beats displayed by acoustic systems. Remarkably, the expression relating the in- 

crease/decrease in oscillatory cycles to the quantum of hysteresis appears to be model 
independent. The cyclic scenario is extended to spatially anisotropic models and it 
is shown that the anisotropy density decreases during successive cycles if ^ pdV is 
negative. 



1. INTRODUCTION 

We live in a universe that is old and very nearly spatially flat. The possibility that these 
two properties of our universe - its age and small spatial curvature - could be related, has 



been the focus of considerable study in cosmology. Tolman wondered whether a progressively 
older universe could be constructed out of repeated cycles of expansion and contraction [l|. 
However he was also well aware of the fact that, for perfect fluids, the equations of motion 
are reversible and so each cyclic epoch is identical to the next. To construct a universe in 
which successive epochs were of longer duration Tolman postulated the presence of a viscous 
fluid. Viscosity leads to an asymmetry in pressure during expansion and contraction which, 
in turn, results in a progression of cyclic epochs of successively longer duration. 

However Tolman did not have a prescription for avoiding the Big Bang singularity and 
had to assume it a priori. An important later development which addressed both the age 
issue and the flatness problem was inflation. By driving the flatness parameter, Q, towards 
unity, inflation ensured that the rapidly expanding universe, even if spatially closed, would 
expand for a long duration of time. However inflation did not address the issue of the big 
bang singularity and it has been shown that although inflation could be eternal in the future, 
its past spacetime is necessarily incomplete [2[. 

The present paper further develops the central idea's of an oscillatory universe and at- 
tempts to synthesise elements of cyclic cosmology with the inflationary paradigm. Extending 
the arguments originally proposed in |3| we demonstrate that a universe filled with a scalar 
field possesses the intriguing property of 'hysteresis". Cosmological hysteresis is related to 
the fact that the pressure of a scalar field is usually asymmetric with respect to expansion 
and contraction: Pcxpansion < -Pcontraction- This asymmetry leads to the development of a hys- 
teresis loop, § PdV 7^ 0, during each oscillatory cycle. The loop can cause consecutive cycles 
to be larger in amplitude and in duration. While the asymmetry between expansion and 
contraction is largest for inflationary potentials, the phenomenon of cosmological hysteresis 
appears to be generic and is observed also in potentials which do not give rise to inflation. 

In § 2 of this paper we develop the equations which relate cosmological hysteresis to the 
amplitude of successive cycles and demonstrate that these equations have a universal form 
which is independent of the scalar field potential responsible for hysteresis. Furthermore, 
the presence of hysteresis appears to be robust, and is shown to exist for quite general mech- 
auisms of siugular.ty avo.dance such as those predicted by Braueworld cosmology |4] aud 
Loop Quantum Gravity J5|, iGj. In §3 we demonstrate that § PdV < for flat potentials, 
which leads to the increase in amplitude of successive cycles. A remarkable feature of this 
scenario is that the universe appears to possess an 'arrow of time' even though the field 



equations are formally time reversible ! For steep potentials the value of ^ PdV can be neg- 
ative as well as positive. In this case the phenomenon of Cosmological hysteresis can adorn 
the universe with quasi-regular oscillations, or heats, resembling those in acoustic systems. 
Section 4 discusses the behaviour of a massive scalar field during cosmological contraction. 
One finds that the field can grow to sufficiently large values during the contracting phase to 
give rise to a long duration inflationary phase at the commencement of the next expansion 
cycle. The behaviour of spatial anisotropy in the cyclic scenario is examined in § 5 and a 
brief discussion of our results is presented in § 6. 

2. COSMOLOGICAL HYSTERESIS 



An oscillatory universe requires two essential ingredients: [A] A mechanism for singularity 
avoidance (when the matter density is high) and [B] a mechanism for inducing contraction 
(when the matter density is low). Below we provide a brief summary of the assumptions 
adopted in this paper regarding both A and B. 



[A] Cosmological Bouncing Scenario's have been widely studied; see 
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17 



18J and 



[19| for a review. Within the framework of general relativity (GR) a necessary condition for 



avoiding the 



by matter 20 



jig bang singularity is the violation of the energy conditions usually satisfied 



21| . An alternative viewpoint considers GR to be an 'effective' theory requir- 



ing modification when the space-time curvature becomes enormous. The initial big bang 
singularity can be successfully replaced with a 'bounce' in theories incorporating both these 



sets of ideas, examples being Branewor 



string theory motivated models [lO, 
action 23l-l26l| . 
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d cosmology j^|22|. Loop-quantum cosmology |5|, |6|], 



16| and other modifications to the Einstein-Hilbert 



While a bounce in the early universe could arise in any one of the above scenario's, 
the central results of this paper will be model independent and so will not depend upon 
the specific mechanism sourcing the bounce. For illustrative purposes we shall consider a 
bounce which is known to arise in Branewor Id cosmology (with a time- like extra dimension) , 
in which the Friedmann equations are modified to 
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Pc) a^ 
(p + 3p)-^(2p + 3p) 

Pc 



(2.1) 



This equation is also valid in Loop-quantum cosmology (LQC) when A; = |5|,|6|; for spatially 
open and closed models the LQC equations are discussed in [Tf. A more general class of 
bouncing equations which accomodates (12. ip and FRW dynamics as special cases is (m > 1) 

r2^ SvrG r^_/p^-l k 



H' = -^P{1 > „ , , o 



d AnG 



(P + 3p) - <{ (3m + l)p + 3(m + l)p W -^ 



, ...... ,. .. . -v- . -.., . , • (2.2) 

a 3 L (^ J \Pc/ 

Eqn (12 ■2p reduces to (12. ip when m = 1, and to the FRW limit when pc ^ oo. In the 
presence of several components contributing to the pressure and density, one needs to replace 
P = EiPi and p = Y.;Vi in (EH) and (Q. 

From (12. 2 p we see that the universe bounces when p = pc, at which point H = and 
a = 4:7rGm{pc + Pc) > (where we neglect the curvature term and assume Pc+Pc > 0). The 
prevailance of the bounce in f l2.ll) & fl2.2p is not linked to the violation of any of the energy 
conditions by matter, but is caused instead by a departure of space-time dynamics from the 
predictions of GR at large values of the matter density.^ One might also note the following 
fairly general phenomenological prescription for singularity avoidance [3| which provides a 
reasonable approximation to (12.10 & (12. 2p 

a— 7- a, a— 7- — a, (j) ^ (f), 0— 7-0. (2-3) 

At small values of the density (p <^ pc) usually associated with late times, higher order 
terms in the density in ( 12. ip & ( 12. 2p can be neglected, and cosmic expansion is described by 
standard FRW equations 

i i 

[B] Cosmological turnaround, at late times, can take place in several distinct ways: 

Bl. If the Universe is spatially closed {k = 1) and the density of matter drops off faster 
than a~^. This has been the conventional approach of making a matter /radiation 



dominated universe turnaround and contract [l|, |27|] . Indeed for a perfect fluid such 



[1] The value of pc is related to fundamental parameters appearing in Braneworld cosmology/Loop quantum 
cosmology. We do not write them explicitely since the precise form of pc will not be required in this paper. 




time 



FIG. 1: The expansion factor of a spatially closed matter dominated universe is described by the 
cycloid (I23]l . 



as dust (p = 0) the expansion factor in a closed universe is time-symmetric and is 
described by the cycloid (see figure [1]) 



a{ri) = A{1 — cos r/) , t = A{ri — sin tj) 



for which the space-time becomes singular at a{t) = 0. 



(2.5) 



B2. The universe will turnaround if, in addition to 'normal' matter with p{t) > 0, one 
postulates a form of matter, p, whose energy density becomes negative at late times. 
Members of this category include: 

(i) p{t) = —A/a^, with A > and n < 2, where n = corresponds to a negative 
cosmological constant |A| = A. (Bl can also be regarded as being a member of this 
category for A = k and n = 2.) 

(ii) Scalar fields with potentials Vi{(j)) oc cos0 and V2(0) = A0^ — m^0^, allow V{(j)) 
to evolve to negative values as the universe expands. These models can therefore 
source cosmic turnaround. Both potentials contain the possibility of giving rise to a 
transiently accelerating universe thereby providing us with interesting candidates for 



dark energy [28| as we 



are discussed in 



29 



30| 



1 as cyclic cosmology. Other models of transient acceleration 
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B3. A novel means of using (12. ip to obtain a cyclic universe was suggested in [31|. These 
authors noted that the density in a phantom dark energy component increases as the 
universe expands, thus making the p^ term in (12. ip relevant both at early and at late 
times. In this scenario, the bounce at small values of the expansion factor is caused 
by normal matter, while the universe turns around and contracts due to phantom DE. 

The bouncing scenario [A] together with either of [B1]-[B3] gives rise to cyclic cosmology. 

In this paper we shall assume that during some period in its history the universe was dom- 
inated by a massive scalar field. The presence of a scalar can make cosmological dynamics 
much more versatile, as we demonstrate below. 

In this case the Lagrangian density and the energy-momentum tensor have the form 

Tij = di(f)dj(j) - QijC , (2.6) 

and, for a homogeneous scalar field, the energy density and pressure are, respectively, 

P=\^' + V{ct^) , P=\^^- y^^) , (2.7) 

and the scalar field equation of motion is 

■■ ■ dV , , 

(P + 3H(P+—r = 0. (2.8) 

The term 3H<p in (12. 8p behaves like friction and damps the motion of the scalar-field when 
the universe expands {H > 0). By contrast, in a contracting {H < 0) universe, 3H(p behaves 
like anti-friction and accelerates the motion of the scalar field. Consequently a scalar field 



with the potential 33| V^ = Vo0^^, A; = 1, 2 displays two asymptotic regimes [3^ 

p ~ —p during expansion {H > 0) , (2.9) 

p c::i p during contraction (iiT < 0) . (2.10) 



In the words of Zeldovich 35j ^^ There is a moral to he learned from these simple calcula- 
tions. The result by and large conforms to Braun and Le Chatelier's principle, which also 
holds in human relations: 

Every system resists outside forces. 



J I l/l/ = p/p 



+ 1-- 



1-- 




COSMOLOGICAL HYSTERESIS 

FIG. 2: An idealized illustration of cosmological hysteresis. The hysteresis loop shown above has 
^ pdV < evaluated over a single expansion-contraction cycle. Rapid oscillations of the scalar 
field will modify this picture, since the equation of state of the scalar field varies between —1 
and +1 during oscillations. Indeed, as discussed later in the text, scalar field oscillations are 
indispensible for providing a non-vanishing value to the hysteresis loop. Oscillations of the scalar 
during expansion scramble its phase making all values of (j) equally likely at turnaround and leading 
to ^expansion / ^contraction and <f pdV / 0. By Contrast only a unique value of (p, namely (/) = at 
turnaround, enables the scalar field to follow its original trajectory in reverse during contraction, 
resulting in ^expansion = ^contraction and <f pdV = 0. Consequently only a very small set of initial 
conditions does not give rise to hysteresis, the latter appearing to be ubiquitous for scalar field 
models which can oscillate. 



The scalar field expands, and a negative pressure (or tension) builds up. If the expansion 
were created by the motion of a piston in a cyclinder containing cf), the tension would de- 
celerate the piston. On the other hand, if the field is being compressed, a positive pressure 
builds up, producing a force opposing the motion of the piston. " 

As Zeldovich suggests, there is nothing really surprising about the behaviour of the scalar 
field in a closed universe. What is surprising, however, is the observation [3] that as the 
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universe contracts, the scalar field does not perform a time reversal and move back up the 
same trajectory down which it descended during expansion. What is seen instead is a lag 
between the trajectories describing expansion and contraction. This behaviour, shown in 
figure [2], is typical of systems displaying hysteresis.^ 

Let us now discuss the effect of hysteresis on cosmological dynamics. We assume, as in 
B2, that the late time behaviour of the universe is governed by the Einstein equation 

H'^Kp-^, (2.11) 

a" 

where k = ^^, and A > 0, n <2, so that the universe turns around and contracts if matter 
satisfies the strong energy condition p + 3p > 0. Note that y4 = A<0, n = corresponds 
to a negative cosmological constant, whereas A = 1, n = 2 describes a universe which is 
spatially closed. 

The scalar field (12. 6 p has no dissipation and therefore provides us with an example of a 
perfect fluid [36] . However unlike other perfect fluids such as dust or radiation, the expansion 
factor for a cyclic universe filled with a scalar field need not display time-symmetric evolution. 
The reason for this is as follows. 

At turnaround the universe stops expanding and begins to contract. Setting if = in 
(EnD we get 

m = ^. (2.12) 

where pt is the density and at the expansion factor at turnaround. The mass ^ associated 
with the volume a^ is M = pa^, therefore at turnaround, nMt = Aa^'^. 

The work done during each contraction-expansion cycle is related to the hysteresis loop, 
f pdV, as follows 

5W=ipdV=l pdV+l pdV. (2.13) 

J J contraction J expansion 

Setting 5W + 5Mt = we get ^ 

- IpdV = 5Mt = -Sal-"" , (2.14) 



[2] Nature is replete with examples of hysteresis. These range from the behaviour of fcrromagnets under the 
influence of an external magnetic field, to control systems, mechanics and even economics ! 

[3] Our results do not depend upon a in M = apa^ and we set a = 1 for simplicity. 

[4] The relationship 6M — —pSV follows from the conservation equation T^.f. — ^ p + 3H{p + p) — 0. One 
might note that formulae (|2.15p & (I2.20p agree with our numerical results for a wide range of parameters. 
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from where we find the following simple expression relating the change in amplitude of 
successive cycles to the value of the hysteresis loop (omax = o-t) 



\3— n 



3— n 



3— n 



5 (a^aj^- ^ ^ aW I -lata'M = -1- d> pdV , (2.15) 



A 



where the hysteresis loop is evaluated over one complete contraction- expansion cycle, namely 

ipdV ■= I '^pcil^ , (2.16) 

ttrnax being the maximum value of the expansion factor in the i^^ cycle; see figure El 



From (I2.15P we find that the change in amplitude of consecutive cycles is sensitive both 
to the value of the hysteresis loop, ^ pdV, and the mechanism responsible for turnaround. 
Two extreme cases correspond to: (i) the negative cosmological constant (n = 0) for which 



Sa: 



3 

max 



-- (p pdV 



(ii) the spatially closed universe (n = 2) for which 



""max — ^max "'max 



K (b pdV . 



(2.17) 



(2.18) 




FIG. 3: The change in amplitude of successive expansion maxima is linked, via (I2.15P & ()2.16p . to 
the hysteresis loop defined by ^ pdV := j^^ pdV . By contrast, the change in amplitude of successive 
expansion m,inim,a is related, via ()2.20p &: (j2.2ip . to the hysteresis loop defined by § pdV := J pdV . 

It is interesting that a companion relationship to (I2.15P can be derived for the change in 
the minimum value of the expansion factor at each successive bounce. Setting M = pa^ as 
earlier, and noting from (12. 2 p that p = pc at the bounce,^ we find 5Mbounce = (^(Pcamin)^- 



[5] We assume that the role of the curvature term can be neglected at the bounce which is a reasonable 
assumption provided matter satisfies the strong energy condition, p + 3p > 0, near the bounce. 
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The work done during an expansion-contraction cycle is 



6W = i pdV = pdV+ pdV . (2.19) 

J J expansion J contraction 

Setting 6W + ^Mbounce = we get 

5(a^in)' = {a^n}' - {aill^y = -^ fpdV , (2.20) 

where the hysteresis loop is evaluated over one compete expansion- contraction cycle, namely 

pdV ■= I "^^^pdV , (2.21) 

'^min being the minimum value of the expansion factor in the i^^ cycle; see figure [31 
Comparing (I2.20p and (I2.15P allows us to draw the following important conclusions: 

• From (I2.20p we find that the change in the minimum value of the expansion factor 
depends upon the value of the hysteresis loop, ^ pdV, and the cosmological matter 
density at the bounce, pc, but is insensitive to the nature of turnaround. It is impor- 
tant to note that an increase in the minimum value of the expansion factor at the 
bounce, determined by (I2.20p . would lead to a corresponding decrease in the curvature 
parameter k/al^^^ at the bounce, making it easier for a spatially closed universe to 
inflate even though the curvature term may have prevented inflation from occuring 
during earlier cycles (see also 3l. l37l|). 

• By contrast, the change in the maximum value of the expansion factor, determined 
by (I2.15p . depends upon the nature of turnaround as well as § pdV, but is insensitive 
to the density at the bounce. In passing one might note that the value of n in (12. lip 
is not restricted to being an integer. Indeed if turnaround is sourced by a dynamical 
dark energy model such as the scalar field with potential V{iIj) oc cos (A^) then, as ijj 
rolls towards the negative minimum of V, kinetic terms will ensure that n never quite 
reaches n = 0, the value suggestive of a negative cosmological constant. Consequently 
one might expect a scalar field induced turnaround to mimic Vo/a" with Vq < and 
-2 < n < 0. 

• Note that in the conventional general relativistic framework described by (12. 4p both 
the bounce as well as turnaround can be sourced by the (positive) curvature term. In 
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this case the condition for the bounce demands that, close to it, matter violates the 
strong energy condition p + 3p > 0. On the other hand the universe turns around if 
matter during much later times satisfies the strong energy condition ! It is easy to 
show that in this case (12.151) is replaced by ^ctmax = ^'^ i^ pdV, with ^ pdV evaluated as 
in fl2.16p . while fl2.20p is replaced by a similar expression 6a^i^ = —k § pdV , but where 
§ pdV is evaluated as in f l2.2ip . Within such a conventional setting bouncing models 
with a scalar field were studied in |32|. However it seems unlikely that hysteresis could 
arise in such a scenario since p ~ — p was found to arise during contraction (inducing 
the bounce) and a similar equation of state arose again during the early stages of 
expansion. So the asymmetry between ^contraction and ^expansion, which is an essential 
requirement for hysteresis, is virtually absent in this case. 

• Note also that in a realistic cosmological model only some of the matter degrees of 
freedom driving cosmic expansion in (12. 2p &; (12. 4p may display hysteresis. In this case 
one might expect (I2.15P and (I2.20p to generalize to 

5 (ainax)^"" = --7 j> PadV , 6 {aminf = j> PadV , (2.22) 

where pa denotes the pressure associated with the matter component displaying hys- 
teresis Zc Po-jexpansion 7^ Pa, contraction- 

3. THE BEATING UNIVERSE 

While the equations describing cosmological hysteresis appear to be model independent, 
the concrete value of the hysteresis loop f pdV is related to the dynamics of the scalar field 
and, in particular, to the form of its potential V{(f)). Let us consider some simple potentials 
which generate hysteresis. 



[1] The 'chaotic' potential V{(f)) = VqcJ) gives rise to inflation for k < few |33|. In 



addition, this potential can also serve to describe 'fuzzy' cold dark matter if fc = 1 and 



Vq is sufficiently small [38|. It is interesting that, depending upon the value of Vq, this 
potential can display a steady increase in the amplitude of successive cycles, as well as other 
interesting features such as heats and stochasticity. Figure H] illustrates how the expansion 
factor grows with each successive cycle when turnaround is sourced by a positive curvature 
term (right panel) and a negative cosmological constant (left panel). The issue of heats and 
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FIG. 4: The presence of hysteresis can increase the amplitude of successive expansion maxima in 
a cyclic universe sourced by the potential V{(j)) = |rn,^0^. In the left panel, cosmic turnaround 



is caused by the presence of a negative cosmological constant, so that n = in (|2.1ip . In the 
right panel, turnaround is caused by a (positive) curvature term, so that n = 2 in (j2.1ip . In both 
cases the increased amplitude of successive cycles is described by the hysteresis equations (j2.15p & 
(j2.20p . The successively increasing expansion cycles appear to endow the universe with an arrow 
of time even though the equations governing cosmological evolution are formally time reversable 
and there is no entropy production. 

stochasticity will be discussed slightly later when we turn our attention to the cosh (A0) 
potential. 

A closed oscillating universe holds interesting consequences for the density parameter 



n-l = {aHY 



(3.1) 



As we just saw, cosmological hysteresis can lead to successively increasing expansion cycles 
which would draw the value of Q towards unity. Indeed, (13. ip clearly demonstrates that, for 
an identical value of H, larger values of a{t) will result in smaller values for fi — 1. Clearly, 
a universe with strong hysteresis {f pdV is large and negative) will require comparatively 
fewer oscillatory cycles to reduce the value oi Q — 1, as compared to one in which hysteresis 
was weaker. Thus in a universe with progressiviely increasing expansion cycles, such as the 
one shown in fig. HI the value of the flatness parameter^will gradually be drawn closer to 
unity gently ameliorating the flatness problem; see also 
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For polynomial potentials V oc (p"^^ the end of inflation is marked by coherent oscillations 
of the scalar fleld and the resulting equation of state is 

(-) - 1^ . (3.2) 

Consequently the density (p^) oc a~^^^^^^^'^ falls off /aster than dust (p oc a~^) for k > 1, for 
instance (p,^) oc a~^ in the case of the A0^ potential. As a result M = p^a^ is no longer a 
conserved quantity for A; > 1 ! Nor for that matter is M = p^a^ conserved during inflation 
(when p0 ~ constant) or cosmological contraction (when p^ oc a^^). The formulae relating 
the quantum of hysteresis to the change in expansion maxima/minima however remain valid 
since they are derived on the basis of the conservation equations, T^.^ = 0, which are 
robust to changes in the form of the inflationary potential and mechanisms for making the 
universe turnaround and bounce. It therefore appears that cosmological hysteresis extends 
considerably beyond the domain of spatially closed FRW models for which it was originally 
explored 3[. 

[2] Another interesting example of cosmological hysteresis and cyclicity is provided by 



the potential {rhp^ = vSttG) 

V = 1/o(cosh Xcp/mp - 1) , (3.3) 



whose suitability for being a dark matter candidate was discussed in 39|, |40 |. 

For X(j) <^1, V (X \^(f)^ and the fleld oscillates as pressureless matter with (p) = 0. 

For A0 ^ 1, V ~ |Vo exp (Xcp/rhp) and the behaviour of the fleld can be assessed using the 

slow roll parameters 

ml (V'V A2 ^2 fy" 



Two extreme cases deserve special mention. 

• If A ^ 1 then {e, 77} <^ 1 and the hysteresis loop § pdV has a large absolute value 
signalling strong hysteresis and a steady growth in amplitude of successive cycles. 

• For A > 1, on the other hand, the slow-roll parameters are large and inflation with its 
associated regimes (12. 9 p & (12.101) need not occur. Surprisingly, | § pdV\ 7^ even in 
this case, and the universe can display striking behaviour for moderate values of the 
control parameter 1 <, A <^ f^w. 



14 



a 




a 




time 



time 





time 



time 



FIG. 5: A cyclic universe sourced by the potential V = Vq (cosh A(/)/mp — 1). Increasing the value 
of the control parameter A causes cyclic behaviour to change dramatically, as illustrated in the four 
panels. In all cases turnaround is caused by a negative cosmological constant: n = in (j2.1ip . 
while the control (steepness) parameter has values A = 2, 2.5 (upper left and right) and A = 4, 6 
(lower left and right). As the control parameter increases the behaviour of the universe undergoes 
a remarkable change. For A <, 2 the amplitude of successive cycles increases (top left), with the 
increase being larger for smaller values of A. For moderate values of A (top right and bottom left 
panels) the universe displays a quasi-periodic pattern reminiscent of beats in sound waves. During 
beats the value of ^ pdV is negative during the first half of the larger (parent) cycle and positive 
during the second half. Beats gradually disappear as the value of A is increased. The universe 
now begins to show oscillatory behaviour in which all cycles have roughly the same amplitude and 
duration signifying ^ pdV ~ (bottom right). In all panels the amplitude of successive cycles is 
governed by the simple formulae (|2.17p and (|2.20p . 
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The expansion factor for the cosh potential is illustrated in figures [5] & |6l Figure [5] 
corresponds to the case when turnaround is sourced by a negative cosmological constant 
while in fig. [6] turnaround is sourced by positive (spatial) curvature. Comparing the two 
figures we see similarities in the behaviour of a{t) as well as interesting differences. In 
both cases, moderately small values of the steepness parameter A give rise to a steady 
increase in amphtude of successive cycles. This arises because the hysteresis loop j pdV is 
a negative quantity which leads to a progressively larger universe. Note that while a spell 
of inflation does ensure f pdV < 0, successively larger cycles do not necessarily imply the 
existence of inflation. Indeed for the scalar field models which we have studied, the presence 
of hysteresis appears to be a rather general phenomenon and all that is required for the 
occurance of successively larger expansion cycles is § pdV < ^ ^expansion < ^contraction- The 
inflationary universe occupies the extreme end of this relationship since ^expansion — —1 and 
^contraction ^ 1, as illustratcd in figured 

As the value of A is increased a fundamental difference can be discerned in the behaviour 
of a(t) in figures [5] and El In figure |5] (second and third panels) one sees a modulation 
in the amplitude of successive cycles suggestive of heats in an acoustic system. Beats in 
the expansion of the universe are characterized by a two fold cyclic pattern with smaller 
duration (daughter) cycles nested within a large (parent) cycle. The origin of heats can be 
traced to periodic changes in the value of the hysteresis loop § pdV. During the first half of 
the parent cycle, § pdV < 0, which leads to a steady increase in the expansion maxima of 
successive (daughter) cycles. As the parent cycle reaches its maximum value the hysteresis 
loop changes sign so that § pdV > during the next half cycle. This leads to a steadily 
diminishing amplitude of daughter cycles during the next (parent) half-cycle in accordance 
with (I2.15p . This behaviour is repeated in a self-similar manner during subsequent parent 
and daughter cycles. As in the case of acoustic beats, the daughter cycle with the smallest 
amplitude is located at the boundary of two parent cycles. But unlike the acoustic case, the 
modulation in the value of a^i^ for daughter cycles is more graded than that in Omax, with 
the latter showing more pronounced changes in amplitude during a given parent cycle. 

The heating universe is robust to small changes in the value of the steepness parame- 
ter. Large values (A >^ few) however lead to diminished hysteresis and an end to the heats 
phenomenon. In this case expansion maxima and minima equalise and the oscillatory be- 
haviour of the universe (fig. [5] lower right panel) begins to resemble that shown in figure [1] 
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FIG. 6: A cyclic universe sourced by the potential V = Vq (cosh A(7!)/mp — 1). Increasing the value 
of the control parameter A causes cyclic behaviour to change dramatically, as illustrated in the four 
panels. In all cases turnaround is caused by a positive curvature term: n = 2 in (I2.1ip . while the 
control (steepness) parameter has values A = 1.5, 2 (upper left and right) and A = 3.3, 6 (lower 
left and right). The top two panels show the expansion factor in logarithmic units. In the top 
left panel each cyclic epoch is larger in amplitude and duration than its predecessor, indicating 
§ pdV < 0. As the steepness parameter is increased the behaviour of successive cycles becomes 
highly stochastic (top right and bottom left) which is indicative of the fact that the hysteresis loop 
can be negative {f pdV < 0) as well as positive {§ pdV > 0) during cycles. Once more we find 
larger amplitude cycles to have a longer duration. Further increase of A leads to a stage when 
hysteresis is virtually absent and all cyclic epochs become similar. In all panels the amplitude of 
successive cycles is governed by the simple formulae (|2.18p and (|2.20p . 



for a cosmological model in which expansion and contraction epochs are identical, so that 

Pexpansion ^contraction &nQ j> puV U. 

The existence of the beats phenomenon appears to depend sensitively on the nature of 
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turnaround. Indeed, it appears that beats are entirely absent if turnaround is sourced by a 
(positive) spatial curvature term. It is interesting that for a closed universe the phenomenon 
of beats is effectively replaced by that of stochasticity, as can readily be seen by comparing 
figure E] with figure [5l During stochasticity larger expansion maxima are interspersed with 
smaller ones. Such a situation arises if both § pdV < as well as § pdV > can occur 
during consecutive cycles, and is the physical basis underlying stochasticity.^ 

The origin of stochasticity lies in the fact that for moderately steep potentials A ~ 
0(1), the pressure of the scalar field during expansion can, on occasion, exceed that during 
contraction. This is especially true when a curvature term is present, since in that case 
H^ = up — k/a'^, and the value of the Hubble parameter is smaller than it would be if 
k/a"^ were replaced by a negative cosmological constant. This situation leads to less friction 
(anti- friction) during expansion (contraction) in the scalar field equation of motion (12.81) . 
altering the nature of hysteresis and accomodating the possibility ^expansion > ^contraction, 
and § pdV > 0. Equation ( 12.18^ then leads to (5amax < 0, in other words the amplitude 
of successive cycles can decrease as well as increase ! This is shown in figure [6] which also 
informs us that the duration of a cycle is related to its amplitude and that larger amplitude 
cycles are of longer duration. (In the universe of figure |5l on the other hand, all cycles are 
of roughly equal duration.) Further increase of the control parameter (A ^ 1) causes the 
potential to steepen to such a degree that expansion and contraction epochs become roughly 
similar, so that ^expansion — ^contraction- In this case hysteresis is virtually absent, § pdV ~ 0, 
and the amplitude of successive cycles equalizes giving (5amax — (figure [6] bottom right). 

Note too that stochasticity is usually present only if the potential is too steep to sustain 
inflation. In the presence of inflation one finds § pdV < 0, with the result that the amplitude 
of successive cycles grows with time and stochasticity disappears. Note too that all of the 
above features, namely: (i) monotonically increasing expansion cycles (for small values of 
the control parameter), (ii) Beats and stochasticity (for moderate values of the control 
parameter) exist also for potentials other than cosh(A0), including V ex 0^^ commonly 
associated with inflation. 



[6] The fact that stochasticity will be absent in the cosmological model shown in figure [5] can be understood 
from the following argument. The equations of motion for this model (in which recollapse is sourced by a 
negative cosmological constant) can be written in terms of cf), cj), H governed by a single constraint. Phase 



space is therefore two dimensional and stochasticity is absent due to topological reasons 



41 
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A complementary means of generating hysteresis and the associated increase in consec- 
utive cycles, is by allowing the fluid filling the universe to be viscous and dissipative [jj. 
If ( is the coefficient of bulk viscosity then the fluid pressure changes from its equilibrium 
value po to p = po — 3(H. Consequently p < po during expansion while p > po during 
contraction. The resulting growth in entropy causes sucessive expansion cycles to be larger. 



The possibility that bulk viscosity might drive cosmic acceleration has been studied in |42| . 
and a recent attempt linking entropy production to increased expansion cycles is discussed 
in 



43|. 



A central difference between this approach and ours is that our system of equations (12. ip , 
(12. 8 p & (12. lip is dissipationless and therefore formally time reversible. Yet the presence of 
cosmological hysteresis {^ pdV ^ 0) endows the universe with a plethora of new features 
including beats, stochasticity as well the possibility of a regular increase in the amplitude of 
consecutive cycles. The reason for this proliferation of possibility rests in the following. 

As pointed out in J3|, the presence of hysteresis is closely linked to the ability of the 
field (p to oscillate. Indeed oscillations appear to be vital for the existence of hysteresis 
since they play the important role of mixing the field in phase-space {(f), 0} due to which 
the value of {<p, (f)} when the universe turns around and contracts is almost uncorrelated 
with its phase space value when the field 0(t) began oscillating. This phase-space mixing 
ensures that (during contraction) the scalar field almost never rolls up V{(f)) along the 
same phase-space trajectory down which it descended (during expansion), thereby ensuring 

-T contraction 7^ -< expansion Slid f pClV J^ 0. 

We therefore conclude that V{(f)) must have a well defined minimum value in some region 
of configuration space in order to allow the possibility of hysteresis. Consider as an alterna- 
tive the potential V{(j)) oc </>"" which does not possess a minimum and is a popular candidate 
for dark energy [4^ . In this case no oscillatory phase is present which will reverse the sign of 
during contraction, causing (pit) to roll up its potential instead of down. As a result V{(f)) 
continuously decreases in value when the universe expands as well as contracts. The role of 
the potential during successive cycles therefore soon becomes negligible, and cosmological 
dynamics becomes solely governed by the kinetic energy 0^ oc a^^. The universe, in this 
case, behaves as if it were dominated by a perfect fluid with the stiff equation of state P = p 
and there is no possibility of hysteresis, as demonstrated in figure [71 
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FIG. 7: The expansion factor of a spatially closed universe dominated by a scalar field with the 
potential V oc (p'~°', a > 0. The universe soon settles into an oscillatory mode in which all cycles 
are of equal amplitude indicating <f pdV = 0. 

4. SCALAR FIELD DYNAMICS DURING CONTRACTION 



An important question which arises in connection with the bouncing scenario is whether 
the scalar field can rise high enough on its potential (during contraction) to provide an 
adequate number of inflationary e-folds during the ensuing round of (post-bounce) expansion. 
This value can be easily estimated as follows. 

As discussed earlier, the inflationary scalar field {V = m^0^/2) in a cyclic universe passes 
through three successive regimes during which: (i) (p ^ mp = G^^^^ and the field slow-rolls 
down its potential resulting in p ~ — p, (ii) (f) <^ mp and the motion of the field becomes 
oscillatory leading to (p) ~ 0. The oscillatory regime exists during expansion as well as 
the early stages of contraction, (iii) During the late stages of contraction the scalar field 
amplitude grows to ~ mp. The field now stops oscillating and begins to climb up its 
potential leading to the stiff equation of state p c:^ p. (This is in response to anti-friction in 
(12. 8 p since H < during contraction.) Below we shall provide a simple analytical estimate 
of how far up its potential the field can climb before beginning its descent, soon after the 
universe has bounced, and the universe starts to expand once more. 

It is useful to note in this connection, the following exact expression describing the motion 
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of a massive scalar field in a spatially flat universe which expands/contracts as a power law 
a{t) oc t^, where p = 2/3(1 + w) and w is the equation of state: 



'^^^^ - ' ^ ^ AJ,{mt) + BY,{mt)\ , u = -^—^ . (4.1) 



mp y a3(t) i "' ^ "' 'j ' 2(1 + w; 

During the oscillatory dust-like phase, u = 1/2, and 

, , , cos (mt + 9) , ^ 

where ^ is a phase constant. Equation (14. 2 p informs us that the amplitude of the scalar 
field, averaged over a certain number of oscillations, increases/decreases in a contract- 
ing/expanding universe as (0^) oc a^^. During contraction, once the field value reaches 
~ mp, the oscillatory regime ceases and the scalar field begins to climb up its potential. 
This marks the commencement of the kinetic regime during which the equations of motion 
begin to be dominated by the kinetic energy of 0, resulting in w ~ 1, z/ = in (14. ip . and 
the solution 

— ^ = AJo(mt) + BYo(mt) , where Jo - 1, 'i^o - - In (mt) , when mt < 1 . (4.3) 
mp IT 

In other words = (pintin^Tit/tin, where 0j„ and tj„ are initial values of the scalar field 
velocity and the cosmic time at the commencement of the kinetic regime^. 

Assuming that the bounce occures at the Planck energy and remembering that ~ m,p 
when oscillations end, we get the following result for the scalar field amphtude at the instant 
of the bounce 

^ "^^ In^. (4.4) 



Vl27r m 
(The actual value is smaller if the bounce occures at energies below the Planck scale.) It is 

worth noting that (14. 4p is valid for any potential which is less steep than the exponential^ and 

in this case the ratio Hp/Hin of Hubble parameters at the bounce and at the commencement 

of the kinetic regime, replaces mp/m in the logarithm of (14. 4p so that 

' "^^In^. (4.5) 



y/Vhi H-u 



[7] This equation can also be derived by noting that <j) + 3iJ0 ~ during the kinetic regipie 
[8] In the derivation we used the massless approximation which fails for steep potentials [45 
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While the above expression is quite general and holds for a fairly wide class of inflationary 
potentials, the number of e-folds during the post-bounce inflationary stage does depend upon 
the concrete form of the inflaton potential. Consider for instance the chaotic inflationary 
potential, V = iv?(fp' /2 with m = 10~^mp, which is in excellent agreement with observations. 
Substituting m = 10~^mp into f l4.4p one obtains (pb ~ 2mp, and an embarassingly small 
value for the number of inflationary e-folds A^ = 27r0^ ~ 25. Surprisingly this does not 
indicate that this model can be ruled out, since, to the scalar field value at the bounce, (pb, 
one needs to add 0^ ^ the 'initial' value of the scalar field at the commencement of the kinetic 
regime. The value of 0j„ depends upon the phase of the oscillatory scalar field and ranges 
from —mp to nip. Moreover, immediately after the bounce the absolute value of the scalar 
field continues to increase despite the fact it now encounters a large amount of friction (the 
post-bounce value of H being positive). Indeed, our numerical simulations^ indicate that 
while the value of the scalar field at the bounce is only 2.5 mp, the field manages to climb to 
the significantly higher value 4.6 nip soon after the bounce, resulting in the observationally 
comfortable value N ~ 130. 

5. THE BEHAVIOUR OF ANISOTROPY IN AN OSCILLATING UNIVERSE 

A central issue in cosmology concerns the class of initial conditions which can give rise to 



a universe resembling our own. This issue was highlighted in 46| which showed that isotropic 



models were a set of zero measure in the space of all homogenous solutions of the Einstein 



equations [47|. Subsequently several mechanisms were identified, including cosmological 



particle creation |48|, which might successfully dissipate a large initial anisotropy within 
short span of time, so that a universe which started out being highly anisotropic would 
rapidly isotropize to a FRW space-time. Perhaps the most successful of these mechanisms 
is inflation, which possesses a no hair proper ty which allows the universe to inflate from a 



fairly general class of initial conditions 49|, |50| . In the context of the present paper we would 



like to ask the question as to how a large amount of anisotropy might impact the behaviour 
of a cyclic universe. For this purpose we shall focus on a Bianchi I universe which expands 

[9] Although our analysis is in broad agreement with that of [10[ our results for N are somewhat smaller than 
those in that paper. We believe this is due to a small typo in [lOJ due to which the value of the scalar 
field at the bounce is overestimated. 
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FIG. 8: The expansion factor (dotted) of a spatially anisotropic Bianchi I universe dominated by a 
massive scalar field. The anisotropy density (solid), given by panis = S/a®, first increases and then 
decreases at the time of each bounce. The hysteresis loop for this universe is negative <f pdV < 0, 
with the result that the amplitude of successive cycles gradually increases. This leads to the waning 
of anisotropy with each successive cycle with the result that Panis becomes vanishingly small at the 
commencement of the fifth cycle and cannot be resolved on the scale of the figure. 

at different rates along the three spatial directions and for which the line element is 

ds'^ = dt^ - Rl{t)dx^ - Rl{t)dy^ - Rl{t)dz^ . (5.1) 

Introducing the directional expansion rate 

H, = R,/Ri, ^ = 1,2,3 (5.2) 

and the mean expansion factor a(t) = (i?ii?2-R3)^ we find the mean expansion rate 

(5^3) 



h='-^ij:h., 



a 3 



j=i 



in terms of which the expression for the anisotropy is simply 

3 v^ 



j:(h„ -Hf=^,. 



(hA) 



a=l 



The resulting (0-0) Einstein equation 



31/2 = SttGp + a + 



(5.5) 
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contains the anisotropy in the RHS as if it were an effective energy density with the stiff 
equation of state Panis = Panis = S/a® [50]. In order to assess the behaviour of anisotropy 
in a cychc scenario we shall assume A < 0, which permits the universe to turnaround and 
contract. As before, we also assume that in the vicinity of the 'Big Bang' extra-dimensional 



effects 



^ modify the FRW equations to 

rr2 SttG ( p\ A E 

3 t Pc) 3 a^ 



I (p + 3p) - ^(2p + 3p) + - - 2- . (5.6) 



a 3 1^ Pc J 3 a° 

Our results shown in figure E] indicate that, if the scalar field is not too heavy {m <^ mpi), 
then the oscillating universe displays monotonically increasing cycles which cause the 
anisotropy density panis = S/a^ to decrease from cycle to cycle. The necessary condition for 
the gradual disappearence of anisotropy and the isotropisation of the universe is therefore 
§ pdV < 0, since this guarantees that successive cycles are larger in amplitude to previous 
ones. Note that this does not require the presence of inflation, although a short burst of 
inflation at the start of every cycle would clearly assist the abatement of anisotropy. 

6. DISCUSSION 

In the present paper we have shown how a cyclic universe containing a self-interacting 
minimally coupled scalar field exhibits interesting relationships, fl2.15p & fl2.20p . linking the 
change in the expansion factor at turnaround/bounce to the net work-done during a given 
expansion-contraction cycle: § pdV . This relationship is quite general and requires only 
that the universe turnaround and bounce, so that cyclicity is maintained. It is legitimate 
to ask whether such a scenario can be realised within the framework of a universe which 
accelerates at late times, and therefore resembles our own. The demand that cyclicity makes 
of cosmic acceleration is that the latter be transient since otherwise our current accelerating 
phase would become a permanent feature of our universe, preventing turnaround. It is well 
known that a transiently accelerating phase can occur in several distinct models of dark 



energy 28l-l3l|. However, an investigation of the presence of hysteresis in such models lies 
somewhat beyond the scope of the present paper and we do not discuss it here. 

Another important issue, left untouched in the present paper, concerns the development 
of density inhomogeneities. The universe today is characterized by structure which is sig- 
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nificantly nonlinear on scales shorter than a few Mpc. The criterion 5p/ p ~ 1 marks the 
onset of nonlinearity and the comoving nonlinear length scale k^]^ can be defined from the 



equality (51 1 



{{6p/p)Y^^ = D+{t)U7i P{k)edkj =1 (6.1) 

where -D+(t) is the growing mode of the linearized density contrast and P{k) is the power 
spectrum of linear density fluctuations. For power law spectra P{k) oc /c" 

k^l^Df'^^'\t), (6.2) 

and, since -D+(t) grows with time, it follows that the nonlinear length scale grows contin- 
uously from small initial values at early times to k^l^ ~ 10 Mpc currently. The growth 
of the nonlinear length scale with time reflects the fact that gravitational clustering takes 
place hierarchically in gravitational instability scenario's based on cold/warm dark mat- 
ter. The power spectrum in such scenario's is characerized by the local spectral index 
Ues = d log P{k) / d log k whose value, in the case of cold dark matter, ranges from rics — 1 
on scales ^ 10 Mpc to n^s — —3 on scales <^ 1 Mpc. From (16. 2 p we find that the growth in 
k^l will be affected both by the change in value of the slope nefr(t) (as successively larges 
scales become nonlinear) and the behaviour of i^_|_(t). The latter is governed, in the absence 
of pressure and on scales significantly smaller than the Hubble length, by the Jeans-type 
equation 

D+ + 2HD+ - 4:7rGpD+ = , (6.3) 

from which we learn that the Hubble parameter behaves like a damping term during expan- 
sion, when H > 0, causing linearized perturbations to grow much slower than they would 
in a static universe^*^: D+ oc exp^/4^TGpt. During contraction, H < 0, the situation is 
reversed since the (negative) Hubble parameter now plays the role of an anti-damping term 
which significantly speeds up gravitational instability relative to expansion. This means 
that the nonlinear length scale will grow rapidly during contraction eventually overtaking 
the comoving Hubble length {aH)~^ which decreases with time in a contracting universe. 
One therefore suspects that a universe resembling ours will become strongly inhomogeneous 



[10] In a realistic scenario which incorporates both DE and turnaround, such as [28], density perturbations 
would be expected to grow as 13+ oc a{t) during matter domination, slowing considerably when the 
universe began to accelerate, and speeding up once more when the universe turned around and contracted. 
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during contraction and it is not clear whether the assumptions which went into the deriva- 
tion of the bouncing equations fl2.ip will hold in this case. (Since most bouncing solutions 



have been derived within the FRW setting and rely on assumptions of homogeneity 19| . 
these concerns are also relevant to them.) 

The above argument necessarily applies to an old universe resembling ours. It need not be 
true for a cyclic epoch of much shorter duration ^^. If a shorter cyclic epoch preceeded ours 
then gravitational instability would have had insufficient time to grow to nonlinear values 



and the universe during contraction might persist in being quasi-homogeneous (see 5^ for 



an analysis of linealized perturbations in cyclic models). Such a situation is illustrated by 



a 



the right panel of figure HI the first three panels of figure [6] and by figure 4 of [3 1 . 

To summarize, we have demonstrated the presence of hysteresis in spatially fiat and 
closed FRW cosmologies filled with a self-interacting scalar field. For these models we have 
established a rather simple relationship between the growth in the expansion factor and the 
quantum of hysteresis, namely (I2.15p . We have shown that, depending upon the value of the 
hysteresis loop, the universe can exhibit a wide format of behaviour including progressively 
larger expansion cycles (when f pdV < 0) as well as stochasticity. An early infiationary 
epoch is not an essential prerequisite for the existence of increasing expansion cycles. What 
is required in this case is that the hysteresis loop be negative, § pdV < 0, and even a small 
asymetry between the fiuid pressure during expansion and contraction can accomplish this. 
In such situations, when ti'expansion < wicontraction but ifexpansion > — 1 we sat thay hysteresis is 
small, in contrast to the situation portrayed in fig [2] when ^expansion — —1 and hysteresis is 
large. As demonstrated in this paper, even in situations with small hysteresis, the moderate 
increase in successive cycles can draw the cosmological density parameter towards unity, 
gently ameliorating the fiatness problem. The same is true for spatial anisotropy whose 
value declines in a universe with growing cycles. The presence of Cosmological hysteresis 
can also adorn the universe with quasi-regular oscillations, or beats, reminiscent to those in 
acoustic systems. One should note here that the phenomenon of beats (and stochasticity) 



[11] The above argument can also be circumvented if, during the bounce, the universe gets 'recycled' so that 
it embarks on its next expansion cycle with a new set of coupling constants and a different inflationary 
potential V {(/)). A related possibility exists within the landscape paradigm if the inflaton potential is 
multi-dimensional since the inflaton field can sample different pieces of the potential during different 



cycles, as pointed out in [11 j 
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is not peculiar to the cosh (A0) potential but has been observed in other potentials as well, 
including the potential V = |r?7,^0^ associated with chaotic inflation. 

Our treatment has focussed on scalar flelds possessing canonical kinetic terms and cou- 
pling minimally to gravity. We flnd it intriguing that while the system which we study is 
fully relativistic, its broad dynamical features can be enapsulated by the well known non- 
relativistic thermodynamic expression 6E = § pdV . It would therefore be interesting to 
explore whether the phenomenon of hysteresis is more general and extends to cosmologies 
in which some of the assumptions of this paper are relaxed, such as the non-canonical scalar 
field models associated with DBI inflation [53|, k-essence 5^, ghost condensate models 55| . 
cosmologies with interacting fluid components ^], etc. One could also ask whether cosmo- 
logical hysteresis exists for scalar flelds which couple non-minimally to gravity, for instance 
through ^R(f)'^, Brans-Dicke and fleld derivative type couplings 57|, or in cosmological mod- 
els featuring non-local gravity [26|. We also leave untouched the interesting issue of quantum 
stability of a cyclic model displaying hysteresis, which touches on issues beyond the scope 
of the present paper 58|]. Finally, it may also be worth enlarging the present analysis to 
(higher-dimensional) anisotropic models in which some of the spatial directions expand while 
others contract, and ask whether a form of hysteresis might exist in this case too. 



Acknowledgments 

We acknowledge useful discussions with Tarun Saini, Yuri Shtanov, Parampreet Singh 
and Sanil Unnikrishnan. 



[1] R. C. Tolman, Relativity, Thermodynamics and Cosmology, Clarendon Press, Oxford (1934). 

[2] A. Borde, A.H. Guth and A. Vilenkin, Phys. Rev. Lett. 90, 151301, 2003. 

[3] N. Kanekar, V. Sahni and Yu. Shtanov, Phys. Rev. D 63 083520 2001 |astro-ph/0101448|. 



[4] Yu. Shtanov and V. Sahni, Phys. Lett. B 557 1 2003 [e-Print: |gr-qc/0208047] . 



[5] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. D 73, 124038, 2006 [e-Print: 
gr-qc/0604013] . 



[6] G. Calcagni, Phys. Rev. D 71, 023511, 2005 [e-Print: gr-qc/041002^; E. J. Copeland, J. E. 



Lidsey and S. Mizuno, Phys. Rev. D 73, 043503, 2006 [e-Print: gr-qc/0510022|; P. Singh, 



27 
Phys. Rev. D 73, 063508, 2006 [e-Print: |gr-qc/0603043 1 . T. Stachowiak, and M. Szydlowski, 



Phys. Lett. B 646, 209, 2007 [e-Print: gr-qc/061012T] ; H. Maeda, V. Sahni, Yu. Shtanov, 
Phys. Rev. D 76, 104028, 2007, Erratum-ibid.D80:089902,2009. [e-Print: 'arXiv:0708.3237j 
[7] A. Ashtekar, T. Pawlowski, P. Singh, K. Vandersloot, Phys. Rev. D 75, 024035,2007 [e- 



Print: |arXiv:gr-qc/0612104| ; P. Singh, F. Vidotto, Phys. Rev. D 83, 064027, 2011 [e-Print: 
larXiv:1012.1307j . 
[8] M. J. Rees, The Observatory 89, 193, 1969. 
[9] A. A. Starobinsky, JETP Lett. 30, 719, 1979 
[10] Y-S. Piao and Y-Z. Zhang, Nucl.Phys.B 725, 265, 2005 [e-Print: |gr-qc/0407027] ; 



[11] Y-S. Piao, Phys. Rev. D 70, 101302, 2004 [e-Print: ,hep-th/0407258| ; Y-S. Piao, Phys. Lett. 

B 677, 1, 2009; J. Zhang, [e-Print iarXiv:1102.5395j : 
[12] P.W. Graham, et fl/.[e-Print{arXiv:1109.0282]. 
[13] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 



2001 farXiv:hep-th/0103239]; J. Khoury, B. A. Ovrut, N. Seiberg, P. J. Steinhardt and 



N. Turok, Phys. Rev. D 65, 086007 2002 |arXiv:hep-th/0 108187]. P. J. Steinhardt and 
N. Turok, Phys. Rev. D 65, 126003 2002 [a rXiv:hep-th/0111098| ; Science 296, 1436 
2002. J. Khoury, P. J. Steinhardt and N. Turok, Phys. Rev. Lett. 91, 161301 2003 
|arXiv:astro-ph/0302012] . J. Khoury, P. J. Steinhardt and N. Turok, Phys. Rev. Lett. 92, 
031302 2004 [arXiv:hep-th/0307132] . 
[14] E. I. Buchbinder, J. Khoury and B. A. Ovrut, Phys. Rev. D 76, 123503 (2007) 
|hep-th/0702154| . 



[15] R. KaUosh, L. Kofman and A. Linde, Phys. Rev. D 64, 123523, 2001 |arXiv:hep-th/0104073] ; 
J. Khoury, B. A. Ovrut, P. J. Steinhardt, N. Turok, |arXiv:hep-th/0105212( R. Kallosh, J. 



U. Kang, A. Linde, V. Mukhanov, JCAP 0804:018,2008 [arXiv:0712.2040 j: A. Linde, V. 



Mukhanov, A. Vikman, JCAP 1002:006,2010 [arXiv:0912.0944J; B. Xue and P. J. Steinhardt, 

Phys. Rev. D 84, 083520 (2011) |arXiv:1106.1416j . 
[16] J. -L. Lehners, Class. Quant. Grav. 28, 204004 (2011) [arXiv:1106.0T72] . 
[17] J. D. Barrow and M. P. Dabrowski, MNRAS 275, 850 (1995). R. Diirrer and J. Laukenmann, 

Class. Quant. Grav. 13, 1069 (1996). 



[18] T. Biswas, e-Print: arXiv:0801.1315 T. Biswas and S. Alexander, Phys. Rev. D 80, 043511, 
2009 [e-Print: larXiv:0812.3T82] 



28 

[19] M. Novello and S.E.Perez Bergliaffa, Phys.Rept. 463 127, 2008. [e-Print: larXiv: 0802. 1634] 
[20] S. W. Hawking and G. F. R. Ellis, The Large Scale Structure of Spacetime, Cambridge Uni- 
versity Press, Cambridge (1973). 
[21] C. Cattoen and M. Visser, Class. Quant. Grav.22, 4913, 2005 [e-Print: |gr-qc/0508045] 
[22] P. Kanti and K. Tamvakis, Phys. Rev. D 68, 024014, 2003 [hep-th/0303'073) ; S. FofFa, 



Phys. Rev. D 68, 043511,2003 [hep-th/030 4004 1 ; J. L. Hovdebo and R. C. Myers, JCAP 



0311:012,2003. hep-th/0308088|; C.P. Burgess, F. Quevedo, R. Rabadan, G. Tasinato and I. 



Zavala, JCAP 0402:008,2004. [h ep-th/0310'T22l ; C. Germani, N.E. Grandi and A. Kehagias, 
Class. Quant. Grav.25, 135004,2008. [ hep-th/0611246] ; D. A. Easson, R. Gregory, G. Tasi- 
nato and I. Zavala, JHEP 0704:026,2007. [hep-th/070r252] ; R. Maier, I. Damiao Soares and 
E.V. Tonini, Phys. Rev. D 79, 023522,2009; 

[23] M. A. Markov, Ann. Phys. (N.Y.) 155, 333 (1984). V. Mukhanov and R. Brandenberger, 
Phys. Rev. Lett. 68, 1969 (1992); R. Brandenberger, V. Mukhanov and A. Sornborger, Phys. 
Rev. D 48, 1629 (1993). 

[24] A. D. Linde, Phys. Lett. B 211, 29 (1988). 

[25] M.R. Setare, Phys. Lett. B 602, 1, 2004 [e-Print: i hep-th/0409055] ; D.H. Coule, Class. Quant. 



Grav.22, R125, 2005 [e-Print: gr-qc/0412026|; F. T. Falciano, M. Lilley, P. Peter, Phys. Rev. 



D 77, 083513, 2008. [e-Print: l arXiv:0802.1196] H-H. Xiong, Y-F. Cai, T. Qiu, Y-S. Piao, 
X. Zhang, Phys. Lett. B 666, 212, 2008. [e-Print: larXiv:0805.04r3 ] M. V. Battisti, Phys. 
Rev. D 79, 083506,2009. Yi Ling, Wei-Jia Li, Jian-Pin Wu, JCAP 0911:016,2009 [e-Print: 
larXiv:0909.4862 ]: X. Zhang, Phys. Lett. B 683, 81,2010 [e-Print: larXiv:0909.4940j : R. Bran- 
denberger, [e-Print: larXiv:0904.2 835] Yi-Fu Cai and E.N. Saridakis, JCAP 0910:020,2009. 
[e-Print: larXiv:0906.r78 9^: Yi Ling and Q. Wu, Phys. Lett. B 687, 103, 2010. [e-Print: 
larXiv:0811.2615] : C. Lin, R. H. Brandenberger and L. P. Levasseur [e-Print: iarXiv:1007.2654j 
Yi-Fu Cai and E.N. Saridakis, Class. Quant. Grav.28, 035010,2011 [e-Print: larXiv:1007.3204] : 
Yi-Fu Cai, R. Brandenberger and X. Zhang, JCAP 1103:003,2011 [e-Print: larXiv:1101.0822] : 
K. Zhang, P. Wu and H. Yu, e-Print: larXiv: 1202. 13971 



[26] S. Nojiri and S. D. Odintsov, Phys.Rept. 505, 59, 2011 [e-Print: iarXiv:1011.0544 j: T. Biswas 
et al, Phys. Rev. Lett. 108, 031101, 2012 [e-Print SXiv:1110.5249 ]: A. S. Koshelev and S. 
Yu. Vernov, e-Print iarXiv:1202.1289] 

[27] Ya. B. Zeldovich and I. D. Novikov, The Structure and Evolution of the Universe, University 



29 

of Chicago Press, Chicago, (1983). 
[28] U. Alam, V. Sahni and A.A. Starobinsky, JCAP 04(2003)002. 
[29] J. Frieman, C.T. Hill, A. Stebbins and I. Waga, Phys. Rev. Lett. 75, 2077, 1995; K. Choi, 

Phys. Rev. D 62, 043509 2000; R. Kallosh, A. Linde, S. Prokushkin and M. Shmakova, Phys. 

Rev. D 66, 123503, 2002. 
[30] J.D. Barrow, R. Bean and J. Magueijo, MNRAS 316, L41, 2000; V. Sahni and Y. Shtanov, 



JCAP 0311 (2003) 014 [e-Print: astro-ph/0202346 1 ; S. Nojiri and S D Odintsov, Phys. Lett. B 



637, 139, 2006 [e-Print: hep-th/0603062] ; A. Shafieloo, V. Sahni and A. A. Starobinsky, Phys. 
Rev. D 80, 101301, 2009, [e-Print: ' arXiv:0903.5T4T ]: F.E.M. Costa and J.S. Alcaniz, Phys. 
Rev. D 81, 043506, 2010, [e-Print: larXiv:0908.4251 J J.C. Fabris, B. Fraga, N. Pinto-Neto 
and W. Zimdahl, JCAP 1004 (2010) 008 [e-Print: larXiv:0910.3246] : A. Shafieloo, V. Sahni 
and A. A. Starobinsky, Annalen Phys. 19, 316, 2010; X. Chen, Y. Gong, E. N. Saridakis 
and Y. Gong, e-Print: larXiv:1111.6743l Yu.L. Bolotin, O.A. Lemets and D.A. Yerokhin, e- 
Print: larXiv:1108.02n3t M. Li, X-D. Li, S. Wang and Y. Wang Commun.Theor.Phys. 56, 
525, 2011 [e-Print: arXiv: 1 103.5870] : M. Campista, B. Santos, J. Santos and J.S. Alcaniz, 



Phys. Lett. B 699, 320, 2011 [e-Print: larXiv: 1012.3943] : A.C.C. Guimaraes and J.A.S. Lima. 
Class.Quant.Grav. 28, 125026, 2011, [e-Print: larXiv: 1005.2986] : 
[31] M. G. Brown, K. Freese and W. H. Kinney, JCAP 0803:002,2008. [e-Print: |astro-ph/0405353| 



[32] A.A. Starobinsky, Sov. Astron. Lett., 4, 82, 1978. D. Page, Class. Quant. Grav. 1, 417, 1984. 

A.Yu. Kamenshchik, I.M. Khalatnikov and A.V. Toporensky, IJMP D6, 673, 1997. N.Cornish 

and E.P.S.Shellard, Phys. Rev. Lett. 81, 3571, 1998. 
[33] A. D. Linde, Phys. Lett. B 129, 177, 1983 
[34] V.A. Belinsky, L.P. Grishchuk, I.M. Khalatnikov and Ya.B. Zeldovich, Sov. Phys. JETP 63, 

195, 1985; Phys. Lett. B 155, 232, 1985; V.A. Belinsky, et a/.Prog. Theor. Phys. 79, 676, 

1988; G. Felder, L.A. Kofman and A. Linde, Phys. Rev. D 66, 023507, 2002. 
[35] Ya. B. Zeldovich, My Universe: Selected Reviews, edited by B. Ya. Zeldovich and M. V. Sazhin, 

Harwood Academic p. 93 (1992). 
[36] V. Faraoni, larXiv: 1201. 14481 

[37] J. E. Lidsey and D. J. Mulryne, Phys. Rev. D 73, 083508, 2006 [e-Print: |hep-th/0601203| ; 
[38] W. Hu, R. Barkana, A. Gruzinov, Phys. Rev. Lett. 85, 1158 (2000) 
[39] V. Sahni and L.M. Wang, Phys. Rev. D 62 103517, 2000. 



30 

[40] T. Matos and L. Arturo Urena-Lopez, Phys. Rev. D 63, 063506, 2001. 

[41] G. Nicolis, Introduction to Nonlinear Science, Cambrige University Press, 1995. 

[42] D.A. Dicus and W.W. Repko, Phys. Rev. D 67, 083520, 2003 [e-Print: |hep-ph/0211109] ; W. 



Zimdahl, D.J. Schwarz, A.B. Balatkin and D. Pavon, Phys. Rev. D 64, 063501, 2001 [e-Print: 



astro-ph/0009353|; S. Nojiri and S. D. Odintsov, Phys. Rev. D 72 023003, 2005 [e-Print: 



hep-th/0505215|; J.C. Fabris, S.V.B. Goncalves and R. de Sa Ribeiro, Gen.Rel.Grav.38:495- 



506,2006 [e-Print: astro-ph/0503362 1 ; I. Brevik and O. Gorbunova, Gen.Rel.Grav. 37, 2039, 



2005 [e-Print: !g r-qc/05040Ql] ; G. J. Mathews, N. Q. Lan and C. Kolda, Phys. Rev. D 78, 

043525, 2008 [e-Print: larXiv:0801.0"85 3]: J-S Gagnon and J. Lesgourgues, JCAP 1109:026,2011 

[e-Print: larXiv: 1107. 1503] 
[43] T. Biswas and A. Mazumdar, Phys. Rev. D 80, 023519, 2009. [e-Print: larXiv:0901.4930] 
[44] B. Ratra, Phys. Rev. D 37, 3406, 1988. 
[45] S. Foster, Scalar Field Cosmological Models With Hard Potential Walls, [e- 

Pr int jgr-qc/9806lT3] . 



[46] Colhns, C.B. and Hawking, S.W. Astroph. J. 180, 317, 1973. 

[47] Note that ambiguities exist in the notion of measure in general relativity as recently 
discussed in J. S. Schiffrin and R. M. Wald, Measure and Probability in Cosmology, e- 



Print larXiv:1202.163 7 , 
[48] Zeldovich, Ya.B. & Starobinsky A.A. (1971) Zh. Eksp. Teor. Fiz., 61, 2161. {Sov. Physics 

JETP, 34 (1972) 1159). 
[49] R.M. Wald, Phys. Rev. D 28, 2118, 1983; A.A. Starobinsky, Sov. Phys. JETP Lett. 37, 66, 

1983. 
[50] I.G. Moss and V. Sahni, Phys. Lett. B 178, 159, 1986. 
[51] V. Sahni and P. Coles, Phys. Rep. 262, 1, 1995. 
[52] D. Wands, Phys. Rev. D 60, 023507, 1999. F. Finelh and R. Brandenberger, Phys. Rev. D 



65, 103522, 2002 [e-Print fhep"-th/0112249] ; J. Martin and P. Peter, Phys. Rev. D 68, 103517, 
2003 [e-Print: hep-th/0307077|; T.J. Battefeld, S.P. Patil, R. Brandenberger, Phys. Rev. 



D D70, 066006, 2004 [e-Print jhep-th/0401010] ; T.J. Battefeld and R. Brandenberger, Phys. 



Rev. D D70, 121302, 2004 [e-Print hep-th/0406180 1 ; T. J. Battefeld and G. Geshnizjani, Phys. 



Rev. D 73, 064013 (2006) |hep-t h/0503160]; T. J. Battefeld and G. Geshnizjani, Phys. Rev. 



D 73, 048501 (2006) hep-th/0506139]; F. C. Mena, D. J. Mulryne and R. Tavakol, Class. 



31 



Quant. Grav.24, 2721, 2007 [e-Print: |gr-qc/0702Q64] ; Y-S. Piao, [e-Print iaFXTv: 1012.2734] : 



S. D. P. Vitenti and N. Pinto-Neto, Phys. Rev. D 85, 023524, 2012 [e-Printl a?Xiv:1111.0888J : 
R. Maier, S. Pereira, N. Pinto-Neto and B. B. Siffert, Phys. Rev. D 85, 023508, 2012 [e- 
Print iarXiv:1111.094 6;: Z-G. Liu and Y-S. Piao, e-Print: aTxrv: 1201. 13711 

[53] G. Dvali and S.H.H. Tye, Phys. Lett. B 450, 72, 1999; J. Garriga and V.F. Mukhanov, 458, 
219, 1999; R. Bean, D.J.H. Chung and G. Geshnizjani, Phys. Rev. D 78, 023517, 2008. 

[54] C. Armendariz-Picon, T. Damour and V.F. Mukhanov, Phys. Lett. B 458, 209, 1999; C. 
Armendariz-Picon, V.F. Mukhanov and P.J. Steinhardt, Phys. Rev. Lett. 85, 4438, 2000; 
Phys. Rev. D 63, 103510, 2001. 

[55] N. Arkani-Hamed, H.C. Cheng, M.A. Luty and S. Mukohyama, JHEP 0405, 074, 2004. 

[56] L. Amendola and S. Tsujikawa, Dark Energy: Theory and Observations, Cambridge University 
Press, 2010. 

[57] C. Germani and A. Kehagias, Phys. Rev. Lett. 105, 011302, 2010; T. Kobayashi, M. Yam- 
aguchi and J. Yokoyama, Prog. Theor. Phys. 126, 511, 2011; C. Deffayet, X. Giao, D. Steer 
and G. Zahariade, Phys. Rev. D 84, 064039, 2011; C. Charmousis, E.J. Copeland, A. Padilla 
and P.M. Saffin, e-Print ia?Xrv:1106.2 000 : S. Tsujikawa, e-Print la?Xrv:1201.5926[ 

[58] A.T. Mithani and A. Vilenkin, [e-Print: arXiv:1110.4096l . 



